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ANALYTICAL  REPORT  ON  DATA  COLLECTED 
DURING  THE  FIRST  THREE  YEARS  OF  THE  COOPERATIVE  PROGRAM 

BETWEEN 

MONTANA  STATE  HIGHWAY  COMMISSION  AND  U.S.  GEOLOGICAL  SURVEY 


Introduction. — In  May  1953,  an  open- file  report  entitled, 
"Floods  in  Eastern  Montana,  Magnitude  and  Frequency",  was  published 
by  the  U.  S.  Geological  Survey  in. cooperation  with  the  Montana  State 
Highway  Commission,  The  report  presented  a  method  of  determining 
the  magnitude  and  frequency  of  expected  floods  applicable  to  any 
area  from  100  to  3,000  square  miles  in  most  of  eastern  Montana,  A 
composite  frequency  curve  was  developed  that  expressed  the  relation 
of  floods  having  a  recurrence  interval  from  1.05  to  20  years  to  the 
mean  annual  flood.  The  composite  curve  was  based  on  the  combined 
frequency  curves  for  16  stations  having  5  or  more  years  of  record. 

An  equation  was  derived  expressing  the  relation  between  mean  annual 
flood,  drainage  area,  and  mean  elevation  of  the  desired  basin. 

The  cooperative  agreement  between  the  Montana  State  High¬ 
way  Commission  and  the  U.  S.  Geological  Survey  stipulated  that  an 
annual  report  be  made  of  the  progress  in  the  collection  of  data  at 
the  crest-stage  gage  stations.  In  addition  to  the  yearly  progress 
report,  it  was  agreed  that  an  analysis  of  the  data  be  made  at  the 
midpoint  of  the  six  year  program.  The  object  of  this  interim  re¬ 
port  is  to  present  the  results  of  that  analysis. 

The  data  collected  have  been  used  to  expand  the  range  of 
the  report  of  May  1953,  mentioned  above,  to  cover  drainage  areas  as 
small  as  one  square  mile. 
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Data  used. — Approximately  40  crest-stage  gages  have  been 
in  operation  since  the  summer  of  1955  for  the  purpose  of  recording 
peak  stages  of  flow.  Of  these,  15  stations  were  considered  for  use 
at  this  time;  the  remainder  were  not  used  because  rating  curves 
have  not  been  defined.  Further  examination  resulted  in  the  elimina¬ 
tion  of  6  more  stations,  mainly  because  it  appeared  that  the  mean 
annual  flood  determined  on  the  basis  of  the  very  short  period  of  re¬ 
cord  was  not  reasonable  due  to  unusual  occurieuices  having  been  ex¬ 
perienced  at  those  stations. 

Records  for  15  stations  used  in  the  original  report  were 
also  used  in  this  analysis.  Little  Beaver  Creek  near  Marmarth,  N. 
Dak.,  did  not  correlate  with  the  remainder  of  the  stations  when  the 
different  basin  characteristics  were  introduced,  and  therefore  had 
to  be  eliminated  in  this  analysis.  A  total  of  24  stations,  15  long¬ 
term  stations  and  9  short-term  stations  were  used  in  the  multiple 
correlation  covering  a  range  in  drainage  area  from  0.B  to  3,000 
square  miles.  All  of  the  stations  used  are  listed  in  table  1. 

Flood  frequency  relations.— Various  techniques  for  deter — 
mining  the  frequency  and  magnitude  of  floods  at  a  single  gaging  sta¬ 
tion  from  past  records  have  been  developed;  however,  only  two  methods 
were  used  in  this  report:  the  partial  duration  series  and  the  annual 
flood  series.  It  has  been  found  that  the  partial  duration  series  and 
the  annual  flood  series  give  essentially  identical  results  for  inter¬ 
vals  greater  than  ten  years.  For  this  reason,  the  annual  flood  series 
was  used  for  the  15  long-term  stations  to  determine  the  mean  annual 
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flood.  The  partial  duration  series  was  necessarily  used  for 
the  crest-stage  gage  stations  to  determine  the  mean  annual 
flood.  The  mean  annual  flood  for  each  of  the  24  stations  is 
tabulated  in  table  1. 

Composite  frequency  curve. — Before  a  group  of  stations 
can  be  combined,  a  test  of  homogeneity  is  considered  necessary  to 
insure  that  all  of  the  records  are  selected  from  a  region  with  uni¬ 
form  flood- frequency  characteristics.  This  test  for  homogeneity 
requires  that  not  only  the  mean  annual  flood  be  determined,  but  also 
the  ten-year  flood  at  each  station  be  determined.  Due  to  the  short¬ 
ness  of  record  at  the  crest-stage  gage  sites,  it  was  impossible,  at 
this  time,  to  determine  the  ten  year  flood  and  therefore  impossible 
to  test  each  station  for  homogeneity.  It  is  felt,  however,  that  at 
the  end  of  six  years  the  ten  year  flood  can  be  determined  from  an 
extension  of  the  flood  frequency  curve.  In  order  to  make  use  of  avail¬ 
able  data  it  was  assumed  that  all  of  the  stations  are  homogeneous. 

The  long-term  stations  were  grouped  together  for  the  purpose 
of  computing  average  flood  ratios  for  each  recurrence  interval.  This 
composite  flood- frequency  curve  is  shown  in  figure  1,  and  is  the  same 
as  that  computed  for  the  original  report. 

A  comparison  between  the  long-term  stations  and  the  short¬ 
term  stations  for  the  common  period  of  1956-58  was  deemed  necessary. 

Of  the  15  long-term  stations,  data  for  only  9  stations  were  avail¬ 
able  at  this  time.  However,  it  was  felt  that  these  were  a  sufficient 
number  of  stations  to  be  representative  for  that  group.  Figure  2 
shows  the  composite  frequency  curves  for  the  9  crest-stage  gage  sta- 
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tions,  the  9  long-term  stations,  and  the  curve  from  figure  1,  all 
plotted  on  basis  of  the  partial  duration  series.  It  can  be  seen 
that  the  long-term  and  short-term  stations  compared  quite  favorably 
for  a  ratio  to  the  mean  annual  flood  above  1,  while  below  1  the 
seemed  to  diverge.  The  divergence  of  the  lower  portion  of  the  curve 
for  the  9  long-term  stations  reflects  a  very  dry  period  in  the  east¬ 
ern  and  northeastern  parts  of  Montana,  Scattered  thunder-storms  often 
cause  some  runoff  at  the  small  area  crest-stage  gage  sites.  These 
same  storms,  however,  may  be  of  such  small  extent  that  their  effect 
on  the  larger  drainage  areas  is  negligible.  Therefore,  it  has  been 
assumed  that  the  composite  frequency  curve  for  the  long-term  stations 
as  shown  in  figure  1  is  applicable  for  all  sizes  of  drainage  areas. 

More  record  at  the  crest-stage  gages  should  prove  whether  or  not  the 
assumption  is  valid. 

Derivation  of  mean  annual  flood. — In  order  to  use  the  region¬ 
al  flood- frequency  curve,  a  study  was  made  of  the  basin  characteristics 
that  might  effect  the  determination  of  the  mean  annual  flood. 

The  drainage  area  was  the  first  characteristic  considered, 
and  was  found  to  correlate  fairly  well  with  the  mean  annual  flood. 
Substantial  error  would  result,  however,  from  estimating  the  mean 
annual  flood  from  drainage  area  along. 

In  the  original  report,  the  mean  altitude  of  the  drainage 
basin  was  found  to  have  correlative  value.  Further  investigation  has 
found  a  more  easily  determined  basin  elevation  figure  which  also  has 
correlative  value.  This  is  the  elevation  at  the  point  seven-tenths 
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of  the  meander  length  of  the  main  stream  course  upstream  from  the 
gage.  The  meander  lengths  of  streams  above  the  gage  sites  were 
found  by  use  of  a  map  measure. 

Another  factor  found  to  have  correlative  value  was  the  in¬ 
clination,  or  fall,  of  the  stream.  The  inclination  of  a  stream  is 
computed  by  subtracting  the  elevation  of  the  stream  at  the  gage  from 
the  elevation  at  seven-tenths  of  the  meander  length  of  the  main  stream 
course  to  obtain  a  difference  in  elevation.  This  difference  is  then 
divided  by  the  total  meander  length  of  the  main  stream  course,  in 
miles,  to  give  a  factor  which  is  called  inclination,  in  feet  per  mile. 

The  direction  of  flow,  or  aspect  of  the  basin,  was  also 
found  to  have  correlative  value.  Different  values  were  selected  for 
different  sectors  of  the  compass,  ranging  from  10  to  L&.  A  drawing 
showing  the  locations  of  the  values  to  be  used  can  be  found  on  the  Kp 
diagram  in  figure  3* 

The  factors  used  in  the  multiple  correlation  for  the  differ¬ 
ent  basin  characteristics  are  tabulated  in  table  1. 

Formula . — The  formula  for  mean  annual  flood  as  derived  from 

the  24  gaging  stations  by  multiple  correlation  is: 

q?  ^3=751. 6A°  *553  e"0-1^  j-G-362  d-0.301 

where  A  =  drainage  area,  in  square  miles 

E  =  elevation  at  7/10  meander  length  of  main  stream 
course,  in  feet 

I  =  inclination,  in  feet  per  mile 
D  =  direction  of  flow,  in  units. 
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The  standard  error  of  estimate  ranges  between  45  percent  and  -31  per¬ 
cent. 

For  ease  of  application,  the  formula  can  be  given  the  gen¬ 
eral  form:  Q2,33  =  KaKeKjKd,  where  the  K  values  are  derived  from  the 
appropriate  K  diagrams  as  found  in  figures  3  and  4.  An  example  of 
the  use  of  the  K  diagrams  is  shown  in  figure  4. 

Use  of  composite  frequency  curve. — After  the  mean  annual 
flood  has  been  computed,  the  magnitude  of  the  flood  for  the  selected 
frequency  is  determined.  Select  the  flood  ratio  that  corresponds  to 
the  desired  recurrence  interval  from  figure  1.  This  ratio  is  multi¬ 
plied  by  the  mean  annual  flood,  and  the  result  is  the  magnitude  of 
the  flood  that  can  be  expected  to  be  equalled  or  exceeded  on  the 
average  of  once  in  the  number  of  years  of  the  selected  recurrence 
interval. 

Limitations — The  results  given  herein  for  computing  the 
magnitude  and  frequency  of  floods  should  be  considered  provisional 
and  subject  to  substantial  error.  On  possible  error  may  be  in  the 
assumption  that  the  composite  frequency  curve  for  the  long-term  sta¬ 
tions  (see  fig.  1)  is  applicable  for  all  sizes  of  drainage  areas. 
Because  of  this  possible  error,  it  is  suggested  that  the  formula  given 
herein  be  used  only  for  drainage  areas  from  1  to  100  square  miles. 

For  drainage  areas  greater  than  100  square  miles,  the  report  entitled, 
"Floods  in  Eastern  Montana,  Magnitude  and  Frequency",  should  be  used. 
As  more  crest-stage  gage  stations  become  rated  and  the  length  of  re¬ 
cord  increased,  these  data  will  be  included  to  add  confidence  to  the 
present  study  and  will  undoubtedly  lead  to  revisions  and  improvement • 
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It  is  further  advised  that  the  elevation  characteristic 


be  limited  to  5 >000  feet  above  mean  sea  level.  Extension  of  any  of 
the  curves  beyond  the  range  of  data  used  in  the  derivation  of  the 
formula  is  not  advised.  Also,  use  of  the  formula  for  streams  out¬ 
side  the  shaded  area  sho-wn  in  figure  1  of  the  original  report  is  not 
recommended. 

The  results  given  herein  should  be  used  only  as  a  guide, 
not  as  an  exact  solution  for  computing  the  magnitude  and  frequency 
of  floods  in  eastern  Montana, 
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Table  lr Factors . used  in . the . multiple,  correlation  of  mean,  annual 

flood  with  basin  characteristics 


Sta 

No. 

•  Station 

Mean 

Annual 

Flood 

<32.33 

Drainage 

Area 

A 

Ele¬ 

vation 

at 

0.7L 

Incli¬ 

nation 

I 

Direc¬ 
tion  of 

Flow 

D 

1 

N.  F.  Musselshell  River  near  Delpine 

115 

38.4 

6,450 

89.2 

10 

2 

Dry  Creek  near  Van  Norman 

5,000 

2,530 

2,840 

7.5 

20 

3 

Rock  Creek  at  international  bdry 

970 

242 

2,970 

14.5 

20 

4 

Wolf  Creek  near  Wolf  Point 

1,320 

244 

2,660 

20.3 

10 

5 

Redwater  Creek  at  Circle 

1,300 

542 

2,670 

8.5 

30 

6 

E.F.  Poplar  River  at  international  bdry 

1,000 

256 

2,610 

5.3 

10 

7 

Poplar  River  near  Poplar 

4,600 

3,070 

2,500 

5.0 

10 

8 

Big  Muddy  Creek  at  Daleview 

1,850 

330 

2,400 

11.2 

10 

9 

Pryor  Creek  near  Billings 

620 

425 

4,360 

25.6 

30 

10 

Soap  Creek  near  St.  Xavier 

370 

111 

3,800 

42.5 

30 

11 

Pass  Creek  near  Wyola 

380 

119 

5,000 

67.1 

30 

12 

L.  Bighorn  River  near  Crow  Agency 

1,750 

1,190 

4,500 

20.8 

30 

13 

L.  Powder  River  near  Broadus 

1,600 

2,000 

3,880 

9.1 

30 

14 

L.  Missouri  River  near  Alzada 

2,050 

780 

3,930 

9.7 

30 

15 

Beaver  Creek  at  Wibaux 

1,400 

351 

3,050 

9.1 

40 

17 

Wheatland  No.  2  near  Harlowton 

62 

2.8 

5,300 

85.7 

10 

19 

Petroleum  No.  2  near  Winnett 

50 

2.32 

2,780 

40.0 

10 

20 

Petroleum  No.  3  near  Winnett 

30 

0.81 

2,780 

60.0 

10 

21 

Disjardin  Coulee  near  Malta 

35 

3.5 

2,500 

62.9 

40 

22 

E.F.  Duck  Creek  near  Brockway 

76 

12 .4 

2,900 

29 .4 

40 

24 

Box  Elder  Creek  near  Plentywood 

52 

9.4 

2,340 

31.1 

30 

27 

Spring  Cr  at  highway  16  nr  Plentywood 

205 

16.9 

2,290 

30.0 

30 

28 

Wets  Creek  near  Billings 

56 

9.14 

4,000 

30.7 

20 

29 

Yellowstone  No.  1  near  Billings 

91 

2.58 

3,800 

15.0 

20 

I 


Xsx/nnB  nsem .  So  .  aoXteXanoo  9XqXXXx/m .  axtt  nX ,  baa#  siocf  asl-X  9Xc*bT 
_ aoxXaxiad-QBisrio  nxeBd  ridxw  booXl _ 


-oaiXCI 
lo  noXX 

woX3 

a 

-iXonl 

noXcten 

I 

-axa 

noXXav 

Xb 

 JY.il 

93BniB*xa 

B9lA 

A 

nB9M 

iBi/nnA 

booX3 

asP 

iioXXbX8  • 

Bi8 

.oM 

OX 

s.es 

OiM.d 

A.se 

exx 

9nXqX90  iB9n  *i9vXH  XXeri8X98euM  .3  .W 

X 

OS 

e.Y 

0A8tS 

oeets 

ooote 

nfinnoM  asV  iB9fl  3*99iO 

s 

OS 

e.Ax 

0Y9,S 

SAS 

OY? 

T£iW  XBnoXXBni9XnX  Xb  3*99iO  3*ooH 

e 

ox 

e.os 

OddtS 

AAS 

ose<x 

XnXo3  IXoW  iB9n  3*99iO  IXoW 

A 

0  £ 

e.s 

OYd,S 

sAe 

ooe<x 

9Xo*iiO  XB  3*99lO  19XBfcfb9fl 

e 

ox 

e.e 

OId,S 

d;s 

000tX 

^ibcf  XBnoXXBmaXnX  Xb  ibvXH  iBlqo3  .3.3 

a 

ox 

0.£ 

ooe<s 

oYote 

OOd.Jl 

ifiXqo3  iB9n  i9vXfl  ifiXqo3 

Y 

ox 

s.xx 

OOAtS 

oee 

oestx 

wbXvbXbQ  Xb  3*99iO  ^bbuM  gX3 

8 

oe 

d.?S 

ode,4i 

esA 

OSd 

BSfllXXia  iB9n  3*99iO 

9 

oe 

e.sA 

ooete 

XXX 

oYe 

19XvbX  .X3  iBen  3*9910  qBo3 

OX 

oe 

I.  Yd 

ooote 

?xx 

ose 

bXo^W  1B9H  3*9910  88B3 

XX 

oe 

8. OS 

ooe,A 

oextx 

oeYtx 

Xoae^k  woiO  iB9n  19VXH  morfsXa  •  J 

sx 

oe 

1.9 

088te 

ooo*s 

OOd.I 

sx/bBOia  iB9n  i9vifl  i9bvro3  .J 

ex 

oe 

9.9 

oe?te 

oev 

oeots 

BhssXA  iB9n  i9vlfl  XiuoeaXM  .J 

AX 

OA 

1.9 

oeote 

xee 

OQA,X 

30JBdiW  Xfi  3*99l0  19VB93 

ex 

OX 

Y.ee 

ooete 

8.S 

sd 

xioXwoXibH  iB9n  S  .oW  brualXfisriW 

YX 

ai 

O.OA 

08Y.S 

se.s 

oe 

XXannXW  iS9n  S  .oH  imreXoictea 

91 

or 

O.Od 

08YtS 

X8.0 

oe 

XXannXW  isan  e  *oW  m/gXoiisa 

OS 

OA 

e.sd 

ooets 

e.e 

ee 

bXXbM  iB9n  99X000  nXbiBtaiG  XS 

OA 

A  .9$ 

oo?ts 

A.SX 

dv 

YBW3*OOia  1  B9n  3*99l0  3*OlKI  .3.3 

SS 

oe 

x.xe 

0Aets 

A.  9 

se 

bOCW^Xn9l3  1B9n  3*99l0  i9bXa  xoa 

AS 

oe 

o.oe 

<*sts 

9.61 

eos 

boo*n{;Xn9Xa  in  61  TCBwrfaXrf  Xb  iO  anXiqS 

YS 

os 

Y.oe 

000 1 A 

A 1.9 

*e 

agnXXXXa  iBan  3*99i0  aXaW 

8S 

os 

o.ex 

ooe%e 

8;.s 

19 

aaniXXia  iBen  I  ,oM  snoiwoXXal 

9S 

o 

o 

CM 


O 

C 


o 

in 


o 


o 

fO 


O 

CM 


O 

<T> 

CO 

ID 

m 


fO 


CM 


in 

5 

rO 

CM 


Recurrence  interval,  in  years 


rar 


’ 


'■  ■  • 

ff‘  +  . 
' 


w-j  .,-i 

yhri 


■ 

■ 

. 


■ 

■ 

■ 


■  • 

n  I m 

/ 


, 


.  .  . 


, 


i  ■ 

■ 


,■ 


' 


'  !: 


' 


4 


I  ~ 


